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Abstract 

We formulate the thermal renormalization group, an implementation of the Wilsonian RG in the real- 
time (CTP) formulation of finite temperature field theory, for fermionic fields. Using a model with 
scalar and fermionic degrees of freedom which should describe the two-flavor chiral phase-transition, we 
discuss the mechanism behind fcrmion decoupling and universality at second order transitions. It turns 
out that an effective mass-like term in the fermion propagator which is due to thermal fluctuations 
and does not break chiral symmetry is necessary for fermion decoupling to work. This situation is 
in contrast to the high-temperature limit, where the dominance of scalar over fermionic degrees of 
freedom is due to the different behavior of the distribution functions. The mass-like contribution is 
the leading thermal effect in the fermionic sector and is missed if a derivative expansion of the fermionic 
propagator is performed. We also discuss results on the phase-transition of the model considered where 
we find good agreement with results from other methods. 
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1 Introduction 



Finite temperature field theory has become a mature part of theoretical particle physics over the past 25 
years. Since the seminal work of Kirzhnits and Linde |]lj] on symmetry restoration at high temperatures 
there has been a strong theoretical interest in methods and results in the field. In recent years, the 
experimental quest for the high temperature state of strongly interacting matter, the quark-gluon plasma, 
has attracted additional interest. Nevertheless, there is a large number of open problems both on a 
phenomenological as well as on a fundamental level. These range from the equation of state of nuclear 
matter at high temperature and density (of interest e.g. in astrophysics) to questions about the approach 
of a quantum system to thermal equilibrium. 

The most general approach to quantum field theory in a nontrivial "environment" (as specified by a 
density matrix at some initial time t{) was proposed by Schwinger 0] and Keldysh || and is known under 
the name of the closed time-path (CTP-)formulation of quantum field theory. This approach allows for 
the calculation of operator expectation values in quite general systems, including systems in thermal 
and chemical equilibrium. The special case of equilibrium physics has, because of its relative simplicity, 
attracted most of the attention in the field. It has also led to two (in principle equivalent) formulations, 
known as the Matsubara- or imaginary-time and the real-time formulation of thermal field theory (see 
e.g. |, §). As the names suggest, in the (more widely used) Matsubara-approach one formulates the 
theory in Euclidean space (with a finite time- interval) , the price to pay for this being the necessity of 
analytically continuing all time-arguments of Green- functions at the end of the calculations. The real- 
time formulation, on the other hand, works in Minkowski space throughout and thus the need to do 
analytical continuations never arises. 

As long as one is mainly interested in static quantities like the free energy etc, there is of course no need 
to consider time-dependence. Since both formulations should yield the same answers for all quantities one 
will usually pick the technically simpler one - which most of the time will be the Matsubara approach. 
However, even in a calculation within a thermal system, time-dependent correlation functions carry 
important physical information to be interpreted within linear response theory Thus for example 
the decay of out-of-equilibrium fluctuations of a given system is governed by the " damping rate" , being 
connected to the imaginary part of the self energy of the corresponding fields. This quantity plays an 
important role in systems approaching the critical temperature of a second order phase-transition. In 
particular it encodes the phenomenon of critical slowing down || - the fact that the equilibration time 
diverges as T — > T c - and governs the initial density of possible topological defects M. It also provides a 
link between fundamental physics and a commonly used phenomenological approach to out-of-equilibrium 



systems using Langevin-equations [||, ||, [LCj with potentially large qualitative impact on the behavior of 



such systems [11]. The fact that realistic systems will be prepared at some temperature and then cooled 
(or heated) "through" a phase-transition during some finite time thus provides a strong motivation to 
study thermal field theory using the real-time formulation. 

As for quantum field theory at vanishing temperature, there are in general no exact analytical methods 
to do calculations in finite temperature field theory. On the other hand, the methods we use in standard 
field theory may be more or less trivially adapted to systems in or close to thermal equilibrium. One 
may in particular tackle a number of nontrivial problems using standard perturbation theory, i.e. the 
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expansion in some small coupling. It is however widely known that due to the emergence of an additional 
dimensionful parameter - namely the temperature T - the usual criteria for the (superficial) convergence 
of perturbation theory are modified. It turns out that even theories which are weakly coupled at T = may 
not be correctly described in perturbation theory close to a second or weakly first order phase-transition. 
The reason is that at a critical point the correlation length of the order parameter field diverges and 
the theory is correlated over all scales. In such a situation the behavior of the perturbative series is 
drastically modified by infrared phenomena. Thus for example perturbation theory predicts a first order 
phase-transition for 0(iV)-symmetric scalar theories in 3 + 1 dimension, whereas we know this transition 
to be of second order fi~2| , [l3|| . Also the damping rate mentioned above in perturbation theory is found 
to diverge at the phase-transition, yielding "critical speeding up" instead of the expected slowing down. 

A method which is well known in statistical physics and has become more widely spread in the particle 
theory community in recent years is the "Wilsonian renormalization group" invented specifically to deal 
with situations involving divergent correlation lengths [14]. The first application of this formalism to 
field theory is due to Polchinksi ([i~5|], see also ]|l6|). The standard way of adapting the method for a 
relativistic field theory works in Euclidean space-time, and there are a number of concrete realizations of 
this approach. These range from the " environmentally friendly" RG ( [jl7| and references therein) via the 
"auxiliary mass method" [18, ^] to the "exact renormalization group" (ERG) of Wetterich These 
methods are straightforwardly adapted to finite temperature calculations in the Matsubara-formalism. 
However, as pointed out above, there are a number of interesting quantities that may not be obtained 
in a direct way from calculations in the imaginary-time formalism and it thus is attractive to devise 
methods for calculations directly in real time that are also able to deal with critical theories. One such 
method is the "thermal renormalization group" (TRG), proposed by D'Attanasio and Pietroni [21] for 
scalar theories and later also generalized to gauge-theories [22]. It was shown in [23] and ^4| that this 
method provides a good tool to study the phase-transition of (at vanishing temperature weakly coupled) 
scalar theories. Also a first calculation of the behavior of the damping rate of a scalar <^ 4 -theory close to 
its critical temperature was performed in [25] and critical slowing down was indeed observed. 

In the present paper we will be concerned with the formulation of the TRG for theories with fermionic 
degrees of freedom. As an example we will study a simple model displaying a spontaneously broken chiral 
symmetry which is restored at high temperatures. This model involves two flavors of chiral fermions and 
an 0(4)-symmetric scalar sector with at T = spontaneously broken symmetry. A Yukawa-coupling of the 
fermions to the order-parameter field then breaks the chiral symmetry in the presence of a nonvanishing 
VEV. This modification of the sigma-model was investigated in |26j and is believed to describe the 
Nf = 2 chiral phase-transition of QCD in the limit where the axial [/(l)-symmetry is completely broken. 
Based on universality arguments the theory should undergo a phase-transition of second order with 
critical exponents of an 0(4)-symmetric scalar model. After the formulation of the TRG we concentrate 
on the universal critical behavior. It turns out that the standard tool in the application of Wilsonian 
renormalization group equations, the derivative expansion, does not work in the fermionic sector, and we 
discuss the implications of this observation for calculations in models with chiral symmetry. A necessary 
ingredient for fermion-decoupling in the framework of the TRG is an effective mass-term which does not 
break chiral symmetry and is obtained from considering "hard thermal loop" -contributions p^ ]. These 
result in a complicated dispersion relation for fermionic fields and we discuss how the propagators may 
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be simplified without loosing the important effects. 

The paper is organized as follows: In the next section we present the formulation of the thermal 
renormalization group equation for theories involving fermionic degrees of freedom. It will be a straight- 
forward extension of the TRG for scalar theories and for details of the practical application including 



necessary approximations we refer the reader to 1 21, 24 1. In section 3 we introduce the model we will 



use. Section 4 contains the discussion of the fermionic dispersion relations and constitutes an important 
part of the present work. In section 5 we then report on the results we find for the universal aspects 
of the critical behavior of the model under consideration. These compare well with results from other 
approaches and prove our method to constitute a good starting point for the nonperturbative exploration 
of dynamical aspects of the phase-transition. Section 6 gives a summary and conclusions. 

2 The thermal renormalization group for fermions 

In this section we discuss the formulation of the thermal renormalization group (TRG) for fermionic 
degrees of freedom. The TRG equation is a Wilsonian renormalization group equation in Minkowski 
space for the dependence of an effective action on an external cutoff scale A, where this cutoff is imposed 



only on the thermal fluctuations [21]. This is possible since in the real-time formulation of thermal field 
theory the propagators contain two parts, one being identical to the (T = 0)-propagator and another 
one being multiplied by the appropriate thermal distribution. We will impose a cutoff on this second 
part. A somewhat more extensive discussion of the closed time-path approach to thermal field theory in 
connection with Wilsonian renormalization group methods may be found in j24j]. We will here assume 
that the reader is familiar with thermal field theory in real time, and only discuss the implementation of 
the cutoff and the corresponding RG-equations. 

As stated above, the suppression of soft contributions to the path integral is implemented by modifying 
the two-point functions. In the operator approach, the propagators are obtained as the expectation 
values of time-path-ordered products of two field operators at different space-time points. Introducing 
for fermionic fields the Fourier decomposition 

fa) = E / (O \ b( - k > S ) e{ ' lkX) + ^ *M fc > S ) eilkX) 

J (2vr)- j 2cj A . L 

WW = E l j2$kr k \ b Hk,s)u\k,s)e^+d(k,s)vHk,s)e^] (1) 



±s 



the propagator depends on the expectation values of products of b, b\ d, and $ . In thermal equilibrium 
at a given temperature T = f3~ l one has 

{b + {k lS )b{k',s'))p = (2w) 3 2u k N(u k )d ss ,6(k-k') 

(d + (k,s)d(k',s'))p = (27r) 3 2Lu k N(Lu k )5 ss/ 6(k-k') 

(b{k,s)b + {k\s')) p = {2Ttf2u k (l-N{u k ))5 ss ,5{k-k') 

(d(k,s)d + (k',s'))p = (2vr) 3 2 Wfe (l - N(u k j\ 5 ss ,5(k - k') (2) 
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where N(u) is the Fermi-Dirac distribution, 

The above expectation values differ from their bosonic counterparts only by the appearance of N(u) 
instead of the Bose-distribution N(u) and by the minus-sign in the last two equations, this being due to 
the anticommuting nature of the fermion fields. The introduction of the cutoff-scale A now proceeds in 
the same way as in the case of scalar fields [21], namely by replacing 

N(io k )^N(io k ;A) (4) 

with 

iV>,;A)^P^ fOT| ! l>>A . (5) 
V ' \ for « A 

This implies that the "hard" modes are in equilibrium at temperature T whereas "soft" thermal fluctu- 
ations are suppressed. We will write 

N(u; k ;A) = @(\k\,A)N(u k ) (6) 

and in practical applications always use a sharp cutoff with 0(|A;|, A) = 6(\k\ — A). The cutoff thermal 
propagator then reads (a and f3 are spinor indices, T c denotes time-ordering on the contour C) 

iS c ,A(x, x)a/3 = (T c ip a (x)^p(x')) p = 6 c (t - t r )S> A (x, x') a p + 9 c (t' - t)S< A (x, x') a(j (7) 

with 

S> A (x,x') a p = {if) a (x)i)p{x')) A p , S< A (x,x') aj 3 = - {^)p{x')il) a {x'))^ . (8) 
Inserting (|l|) and @) with the modification @ thus gives 

/d^k 
27r5(k 2 - ml){$ + m ) a(3 e- lk ^- x > 

9 c (t — t')9(ko) + c (t' — t)9(—ko) - JV(|feo|)0(|fc|, A)j . (9) 

Now we may adopt the matrix-notation which has become common in real-time thermal field theory (see 
II) and write the fermion propagator as 

with 

A f = $ -mo + ie]" 1 , Af = [ft - m - ie]' 1 . (11) 
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where we have already switched to a ^-function cutoff. Together with an analogous expression for the 



modified scalar propagator [21 we then define the scale dependent partition function Z\[J,rj,r]] as 



Za[J,v,v] 



8<p8ip8ip exp<( i[ -cp ■ D K l ■ <p + ip ■ S K l ■ ip + 



+ J ■ (f> + r] -ip + ip -t] + S iTl t[(p,ip,ip] 



(12) 



From the way the cutoff was introduced, it is obvious that 

Z A ^ 00 [J, V ,fj] = Z^[J,rj,fj] 



(13) 



and that in the opposite limit A — > the usual generating functional for the theory at temperature T is 
obtained. Finally, since all scale-dependence of Z A is explicit and contained in the terms bilinear in the 
fields, it is straightforward to obtain an expression for the A-dependence of Z\. 

As in the Euclidean formulation it is however more convenient to work with a slightly modified version 
of the effective action, the generating functional for the lPI-Green functions. This is obtained from the 
Legendre transform of the logarithm of Z A according to 

1 



-HtlZ\ — J ■ If — Tf ■ t/) — ip ■ rf — -<p ■ D A ■ (f — Ip ■ S a - ip 



Here the fields are now the "classical" fields as (p = s (~ 1 ^ Za ^ e ^ c The sources on the right hand side are 
to be regarded as functionals of the fields, given by the inversion of the above expressions, and are thus 
scale-dependent. The RG-equation governing the scale-dependence of T\ is now easily obtained and one 
finds 

-l 



(14) 



Ad A r A [cD] = ^STt^A^a^ 1 ) 



+ V7 1 



which is a straightforward extension of the scalar TRG-equations | 

$ = ((p,ip,ip) ; $ = ((p,-ip,ip) 



$ denotes a superfield 



(15) 



(16) 



and the supertrace takes into account the anticommuting nature of the Grassmann- fields. Note that we 
have suppressed all indices: The fields carry a spacetime-index, a thermal index, possible internal indices 
and (for the fermions) Dirac indices. All these are traced in (|l5|). Also by our conventions the bare 
inverse propagator-matrix reads 



V A \k) 



( D A \k) 



Sl\k) 



\ 



{Sl\k)) T ) 



(17) 



and we have 



( * 2 r A 



StpScp 
\ StpScp 



s 2 r A 

S<p8ip 

_£r A 

SipSip 
S 2 T A 
SipSip 



5!£a \ 



Sip Sip 
S 2 V A 
SipStp / 



(18) 
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The general properties of the TRG-equation are extensively discussed in [21, 24 1 and we will here only 
make a few remarks in order to make the present work self-contained. Let us first note that no approxi- 
mations have been made at this point. The TRG-equation ( |l5| ) is an exact relation on the same level as 
for example the Schwinger-Dyson equations. However, it is of course impossible to give a general solution 
to this functional differential equation - note that the matrix ( + P7 1 ) is the full propagator. The 



equation is formally a one-loop equation and we may introduce the derivative d\ which is defined to be 
a derivative acting only on the explicit scale-dependence in the bare propagator. Using this derivative 
we have 



Ad A r A [<D] = -STrAS A In [T%> + T>X ) ( 19 ) 

and the one-loop character is apparent. We may thus give a simple recipe for calculations in the framework 
of the thermal renormalization group: To derive the flow-equation of a specific coupling, write down the 
one-loop diagrams that would contribute to this coupling, taking into account all vertices allowed by 
the symmetries of the theory under consideration. The flow-equation is then obtained by considering 
only those parts of the diagram which involve thermal distribution functions (since the A-dependence 
only comes with those in the TRG-approach) , replacing the distribution functions according to (Q) and 
similarly for the Bose-distribution, taking the derivative with respect to the explicit scale-dependence 
and then performing the loop-integral. It is obvious that this integral will be finite. The UV-behavior 
of the thermal parts of any Feynman-diagram is never problematic since any UV-divergences of the full 
theory at finite temperature reside in the (T = 0)-parts of the diagrams. The behavior in the infrared, 
which might be problematic in perturbation theory, is regulated by the introduction of the external scale 
(see (pi)). The formal exactness of this approach is due to the use of full propagators and vertices. 



In order to simplify the calculations somewhat we may make use of an observation made in [28], 
namely the fact that one may get rid of the fields on the lower part of the CTP-contour (the "thermal 
ghost fields") if one is interested in physical observables. To this end one may introduce a functional F 
which depends on one (super) field $ (without thermal indices) and is defined through 

f[*] = r[*i,* 2 [*i]]| Sl=s (20) 

Here $2 is the solution to its equation of motion for arbitrary $1, i.e. we have 

ar[*i,* 2 ] 



= (21) 

$2=<M<E>i] 



For space-time independent fields one has ^f'&i] = $1 and (^) is up to a constant equivalent to 

<$f[*] <&r[*i,* 2 



5$ <5<&i 



(22) 

$ 1= $ 2 =$=const. 



This is the definition that we will use below. The TRG-equation for F^ (which is obtained from ( |2"2| ) 



after replacing F with Ta) is readily found from ( |lq ) and (|22|), Some care is required in the presence 
of fermionic fields, since derivatives with respect to fermionic fields do not commute with a supermatrix 
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like the full propagator. We have for example 

_S_ r (2) ^ r A } | r (2) 8 

Sep A Sip A Sip 

£if> - f + «rfM^ (23, 

where ip (if)) is a bosonic (fermionic) field and M = diag(l, —1, —1) (in the present case). The supertrace 

5r (2) 

is cyclic for supermatrices. Note also that the fermionic derivative of a supermatrix (e.g. ^ ) is not a 
supermatrix. 

In the following section we will first introduce the model that we will be interested in and then turn 
to a discussion of the approximations we have to make in order to convert the functional differential 
equation (|l5|) into a solvable system of (partial) differential equations. 

3 The chiral quark-meson model 

We will be interested in a simple model involving only scalar and fermionic fields with scalar self-coupling 
as well as Yukawa-couplings. It is well accepted that the chiral phase-transition of QCD in the limit of 
two massless flavors and an explicitly broken axial U(l) should be of second order, with 0(4)-model 
universal behavior. This could in principle be modeled by a simple 0(4)-symmetric scalar theory with 
(at T = 0) spontaneously broken symmetry. Recently, Jungnickel and Wetterich have pointed out that a 
modification of the sigma-model might actually be quantitatively connected to QCD by a combination 
of analytical and numerical calculations [26]. This model, the so called "chiral quark- meson model", 



includes (in the limit of strongly broken U(1)a) the usual sigma-model fields, which are however coupled 
to N c dublets of fermionic fields via a Yukawa-coupling. As was pointed out above, perturbation theory is 
not applicable even for small (T = 0)-couplings at this phase-transition. On the other hand, the Wilsonian 
RG was introduced to deal with such situations. The chiral quark-meson model is thus, apart from being 
interesting phenomenologically, an ideal testing ground for our method. 

We shall not discuss in detail here how the model may be connected to QCD and refer the reader to 
|26[ | for details. We will actually in the present work not even try to make quantitative contact with this 
effective description by choosing the values of the (T = 0)-couplings accordingly - this could in principle 
be done and we briefly comment on this below. For the main part of this paper we shall discuss the 
universal physics obtained for this model in different approximations. These are independent of the initial 
values of the couplings and thus also apply to the Nf = 2 chiral phase-transition. Nonuniversal quantities 
will only be considered for small (T = 0)-coup lings in the present work. 

The initial condition that we will use for the effective action T\ as A — ► oo, i.e. the effective action 



of the (T = 0)-theory is (remember that the free two-point functions are subtracted in T\, c.f. (14)) 

(24) 



Ia^oo = d x 



9(T=0) 2 I 7 / V + ^ , ¥ - ^\ , 

It involves a fermion flavor dublet as well as a scalar multiplet containing four real degrees of freedom, 

V> = (V>\V> 2 ) , V = l -(a + mf) (25) 
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Furthermore we use the 0(4)-symmetry to write 

p = Tr(<p^) = ^(a 2 + n 2 ) (26) 

The fermions come with an additional color index and closed fermion-loops correspondingly yield a factor 
of N c . 

The effective action given in ( p4|) is of course only a crude approximation to a full effective action of 
a (T = 0)-theory - one could e.g. include perturbative contributions. We will nevertheless use this simple 
form and report nonuniversal results only for small values of the initial couplings grr=Q) an d h(T =0 \, 
such that the higher order contributions should be small. The reason for doing so is that the universal 
behavior is unaffected by these approximations and that on the other hand the nonuniversal quantities 
would be most interesting using the (T = 0)-parameters that correspond to the realistic case where 
observables calculated from the model (24) match those found in nature. A discussion of this matching 



and the resulting values for the (T = 0)-couplings as obtained from Euclidean versions of Wilsonian 
renormalization group equations may be found in |2(|. One finds for example that the Yukawa-coupling 
is h? ~ O(10) and thus by no means small. If we were interested in nonuniversal features of the chiral 
phase-transition as obtained from the quark-meson model we would have to address the question of initial 
values carefully. Studies using the Matsubara-formulation of thermal field theory in connection with the 



average effective action yield reasonable predictions for e.g. the critical temperature [29]. 

Within our approximations the effective action ( p4[ ) thus constitutes the boundary value for the 
solution of the TRG equation Qi~5D (combined with the definition of Ta, (p0|)). We now have to discuss 
the approximations we will use in order to bring this equation into a solvable form. The simplest 
approximation one may think of is to assume that the effective action at finite temperature always has 
the same form as given in ( p4| ) and that only the couplings - in this case the scalar massterm, the 
quartic coupling and the Yukawa-coupling - depend on the scale A. One then obtains the flow-equations 



of those couplings by plugging T\ from (24) into (15) and matching coefficients on both sides. This 



simple approximation will give completely misleading results, for reasons to be discussed in the next 
section. In the absence of fermions (for h = 0) this approximation would give a second order phase- 
transition, however the values for the universal critical exponents will differ strongly from known results. 
The situation can be considerably improved if one allows for a general form of the effective action for a 
constant scalar field p(x) = p |23|, The effective action for constant field is connected to the effective 
potential through 



f A [<p = y/fy, V = 4> = 0] = (m 2 p - U(p)) J d A x (27) 

for a constant configuration p with the bare mass m. The effective potential corresponds to the free 
energy for constant fields and it approaches a universal form at a second order phase-transition. For the 
results for this universal form in scalar 0(iV)-models, using the approximation of a standard kinetic term 
(as in (pi|)) and a general effective potential, see |^4j| . The critical exponents one finds in this case are 
in much better agreement with the known values than in the case of a simple polynomial approximation 



for U [21]. 



Neglecting corrections to the kinetic term may be motivated by arguing in the framework of an 
expansion of the effective action in powers of derivatives. In this expansion, the effective potential is 
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obviously the lowest order term. The leading term beyond this approximation then is obtained from the 
scalar wave-function renormalization, i.e. using the approximation 



d 4 x 



(Z A - l)Tr (dv^dpip 



+ m 2 p ■ 



Uk(p) 



+ Fermion — contributions 



(28) 



Here the wave- function renormalization Z A , which could in principle be a function of the field, is evaluated 
at the minimum of the effective potential. This will be the approximation we use in the scalar sector 
of the model. It is widely used also in Euclidean versions of the Wilsonian renormalization group, see 
e.g. psjl - In particular we will neglect all imaginary contributions to the two-point functions. These 
would only indirectly affect the scale dependence of the effective potential through a contribution to the 
anomalous dimension. Since the anomalous dimension will turn out to be very small in all situations 
that we will be interested in, we believe that dropping imaginary contributions results in negligible errors 
in quantities like the critical temperature and the critical exponents (with the possible exception of 77, 
see below). On the other hand the imaginary parts are interesting objects by themselves and we plan to 
study them in future work. 

Let us now turn to the sector involving fermionic fields. Here the situation is more complicated. The 
reason for the complication is that in contrast to the bosonic sector, where the leading finite temperature 
effects show up in a local contribution to the mass oc T in lowest order derivative expansion, there is no 
local massterm generated for the fermions. We will come back to this issue in section 4 and for now only 
consider the interaction term. Here one obtains to lowest order in the derivative expansion a general 
function of the fields (which should respect the SU(2) x SU(2) x [/(l)-symmetry of (p4|)). However, one 
expects the fermions to decouple from the infrared physics and universal quantities should be independent 
of the approximations in the fermionic sector. The dependence of nonuniversal quantities as the critical 
temperature should in principle be calculable in perturbation theory and we may, as long as we restrict 
ourselves to small values of hrr=o) an d 9(t=o) > always argue that any induced coupling involving fermions 
are of higher order in the (T = 0)-couplings and thus negligible compared to the Yukawa-coupling. We 
will thus approximate the interaction between fermions and bosons by the original form as given in ( |2"4| ) 
for all values of A. In particular, we always evaluate the Yukawa-coupling at p = and at vanishing 
energies and momenta. 

Let us then discuss as an example the flow-equation for the p-derivative of the effective potential. 
This is obtained from the flow-equation for the (/^-derivative of T A at constant cp and vanishing xjj and tp: 



Ad A {pU'{p)) 



Ad A 



dp 



K, A 



i^=const. ,ip=ip=0 



5p5$5<$> 



(29) 



</?=const. ,ip=ip=0 



Here we have introduced the kernel 



/Ca 



«5 2 r A 

5$5<S> 



+ V 



-1 



-1 



(Ad A v- A 



( 6*T A 
\6<j>6<f> 



-1 



(30) 
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The scalar kernel may be rewritten as [21] 



-K8{\k\ — A)Aa(A* 



= Ag-(A 



b\* 

o) 



N(\h 




(31) 



The fermion-kernel is obtained from ( |3l|) by replacing the bare and full bosonic propagator Ag and A^ by 
the corresponding fermionic ones and substituting the Bose-distribution N(lu) by the minus Fermi-Dirac 
one. Note that we have already inserted a sharp cutoff for (|3l|). Due to the thermal matrix-structure of 
the kernel and the fact that l2£ 



S 3 f A 



r A 



(32) 



<3>2 =< 3?i = < E > 



the trace over the thermal indices in (29) allows us to rewrite the right hand side in terms of Ta- 
Furthermore for real self-energies we may write 



A A (A A r A;; 



Ag(A§)' 



2mp b A (q) 



(33) 



with the spectral density p A and similarly for the fermions. The spectral densities may be obtained from 
the ansatz for the effective action F. In our case, using the approximations discussed above for the scalar 
fields and for the moment dropping any momentum-dependent contribution other than the standard 
kinetic term for the fermions, we find 



2mpUq) 



1 



1 



Z A q 2 -M 2 +ie 

A + m xS B 
V 2 ~ m \sB + ie 



Z A q 2 -M 2 -ie 
A + m x SB 



in 



xSB 



(34) 



For e — ► 0, the spectral densities 



where m^SB = h 2 p/2 and M 2 are the eigenvalues of the matrix ® 2 q . 
yield ^-functions which allow us to perform the go-integration implicit in the trace in (29). The angular 
integral is trivial for constant fields. Finally from the A-derivative of the cutoff, we have another 5-function 
(c.f. (|3l|)) which allows us to perform the integration with respect to the modulus of the loop-momentum. 
The result for the flow-equation of the effective potential is 



Ad A U'(p) 
with 



4tt 2 



3U"( P ) + 2 P u^\p)) ^e(^) + 3U"( P )^e(ol) 



N, 



A 3 



7T 



^35) 



Z A A 2 + U'{p) + 2 P U"{p) 

Z A A 2 + U'(p) 

. 2 h 2 p 
A + —f- 



(36) 



The 0-functions arise from the bosonic spectral density [21] and will be understood in the following. The 
structure of (E35|) is again easy to understand from the one-loop character of the flow-equation: the first 
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contribution is from the radial mode which is massive at the minimum of U(p), the second contribution 
is from the 3 Goldstone-bosons which are massless at the minimum, and the last contribution is from the 
2N C fermionic fields which are degenerate. 

Similar flow-equations may be obtained for the Yukawa-coupling and the scalar wave-function renor- 
malization constant. Before we come to an investigation of the critical behavior of the model using these 
flow-equations in section 5, we will in the next section turn to the question of approximations to the 
fermionic two-point function. 



4 Fermion decoupling, universality and hard thermal loops 

In this section we will discuss how the fermions decouple from the critical behavior of the quark-meson 
model, leading to a second order phase-transition with 0(4)-model critical behavior. We will argue that 
the mechanism of fermion decoupling at a second order phase-transition is different from the mechanism 
behind fermion-decoupling at high temperatures. 

Let us start with discussing the high-temperature limit of our model. We will work for the moment 
with the simplest approximation to the full flow equation for the effective action and only take into 
account the effective potential and the Yukawa-coupling and set Z A = 1. The flow-equation for the 
Yukawa-coupling reads within these approximations 



Ad A h 



A 3 - s 1 
4vr 2 U'(0) 



N(Cu ) N(A) 



A 



where 



o 



A 2 + U'(0) 



(37) 



(38) 



Considering now the limit of T large compared to all scales in the theory we may use the asymptotic 
expressions 

T 



N(x) 



(39) 



and 



N(x) 



(40) 



to obtain the system of flow-equations ( |3~5| ) and (37) in the high-temperature limit 



Ad A h\ 



Ihigh-T 



TA 3 
h 6 



3U"(p) + 2pU^(p) U"(p) 
A 2 + U'{p) + 2pU"(p) A 2 + U'(p) 



+ 



1 



4^ 2 



+ 



(41) 



U'(0) (A 2 + U'(0)) 

The contribution of the fermions to the flow of the effective potential thus vanishes in the high-temperature 
limit, due to the different behavior of the distribution functions (|39| ) and (|4fj|). This is what is referred 
to as "fermion decoupling" and the TRG thus reproduces the fact that the dynamics of the soft modes 
in the high-temperature limit is independent of the fermionic content of the theory. 
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Let us now turn to the question about fermion-decoupling at a second order phase-transition. One 
is interested in this case in a (three-dimensional) scaling solution for the free energy (i.e. the effective 
potential). In order to investigate this scaling solution it is useful to introduce rescaled variables according 



K 



p 

AT 



u(k) 



u(p) 

A 3 T 



(n) j^n—3rpn—ljj(n) 



(42) 



At this point we do not rescale the Yukawa-coupling, since if we assume the fermions to decouple there 
is no reason why this coupling should approach a scaling solution. Introducing furthermore the dimen- 
sionless flow-parameter 



A 



A 
T 



we obtain for the flow of u'{k) from (|3 

\d\u'{n) = -2u'(k) + ku"{k) - 
with the dimensionless quantities 



A 
4^2 



3u + 2ku 



LOrj 



+ 3u 



// n(Aa; 7r ) 



UJ-n 



1 



(43) 



(44) 



1 + u'(k) + 2ku"(k) 
1 + u'(k) 

h 2 K 



1 + 



2A 



(45) 



and n(x) = (exp(x) — 1) 1 , h{x) = (exp(x) + 1) . The flow-equation for the Yukawa-coupling is rewrit- 
ten accordingly as 



Xd\h = h c 



1 1 

4^2^(0) 



n(Acj ) 



+ n(l) 



(46) 



where ujq = 1 + u'(0). We will see that the problems with decoupling within the present approximation 
concern the region of small k, thus to simplify the discussion and to make our point clear we expand the 
effective potential in powers of k. Introducing couplings according to 



n n 



nl 



(47) 



we obtain flow-equations for the u n from (flip. Let us for example consider the scale dependence of ui, 
which is directly obtained from (Q) by setting k = 0. It involves the scalar massterm u\ itself as well as 

2 _ 

the scalar four-point coupling U2 (remember p = ^-) and the Yukawa-coupling h. A scaling solution for 
the scalar couplings u n may arise only in the three-dimensional limit, i.e. for A — > 0. In order to identify 
the corresponding fixed-point values of the couplings we thus expand the flow-equations. In the case of 
«i we have 

3 u 2 „ 1 



\d\U\ = —2u\ 



2vr 2 1 + 



N, 



c n 2 



h z + 0(A) 



(48) 
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and the fermionic contributions obviously do not decouple - a possible fixed-point value would depend 
on h and iV c , in contradiction with the assumption of universality! It is not hard to see what spoils 
the high-temperature argumentation in the present case: It is the fact that we should really not com- 
pare the behavior of the distribution-functions alone, but rather the contribution of couplings times the 
distribution-functions as A — > 0. Then it is clear that for a scaling solution (i.e. at a second order phase- 
transition) while the scalar distribution function diverges for small A as A -1 , the four scalar-coupling 
multiplying it is renormalized to with the same power of A. The product remains finite for A — > 0. In 
contrast, the fermion distribution function does not diverge, but on the other hand the Yukawa-coupling 
is not driven to with A - the product also remains finite and is thus for all A of the same order as the 
scalar contribution. 

This is a particular feature of the critical point. For T ^ T c , away from the second order phase- 
transition, the scalar coupling is not renormalized to zero and indeed the product of coupling times 
distribution function for large — dominates the fermionic contribution as shown above. 

Also notice that the presence of a fermion mass-term (as for k / 0) changes the situation. While 
the behavior of coupling times distribution function does not change (n still approaches a constant), 
now the factor lo^ 1 in the fermion-contribution to \d\u'(K) behaves as J and suppresses the fermion- 
contribution. This observation will help us to see how universality is saved and fermion-decoupling works, 
a question to which we turn now. 



To this end, let us turn again to the truncation used to derive (35) and (pTj). We have pointed out 
that we have neglected all terms beyond the effective potential and the Yukawa-coupling. In a purely 
scalar theory, it is known that using only the effective potential one already obtains quite reasonable 
results for the critical behavior (see e.g. |24|]). This is due to the fact that for the critical theory one is 
effectively using an expansion in the (small) anomalous dimension. There is however one other important 
aspect of the thermal effects on scalar fields which is fully covered by the approximation with only an 
effective potential - the generation of a thermal mass. One may easily check that within the present 
approximation one reproduces the "hard thermal loop" (HTL-) result for the scalar mass term: Consider 
the present model for simplicity with N c = 0. The flow-equation for the mass then reads from above 



Ad A U'( P ) = 



3U"(p) + 2pU^(p)) + W"(p) N{ ^ ] 



(49) 



Adopting a perturbative point of view, neglecting U^\p) and replacing in the spirit of the HTL- 
approximation uj a = = A, we obtain 

dAU'(p) = -6A^ (p) _J__ . (50) 

Remembering that the boundary condition for the flow is the (T = 0)-effective action and neglecting 
- again to leading order in perturbation theory - the scale-dependence of U"(p) we easily solve the 
flow-equation (|50|) to find the temperature-dependent massterm to be given by 

U'(p, T) = U'(p, T = 0) + E^It 2 (51) 

reproducing the leading order perturbative result. The positive contribution cx T 2 of course is also the 
reason for the restoration of the spontaneously broken symmetry. The reason why this leading order 
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p+k 

Figure 1: The one-loop contribution to the fermion self-energy. 



high-temperature result is reproduced in the simple truncation with only the effective potential is that 
it is obtained as a local contribution to the effective action. Of course, this local contribution is fully 
resummed if we solve the flow-equation (^) without further approximations. 

In contrast to this simple behavior of the leading thermal contributions in the scalar case, for the 
fermionic fields there is no local massterm generated by hard thermal loops. It is however well known that 
the leading order corrections to the fermion propagator lead to modified dispersion relations which for 
small momenta resemble the ones of a massive particle without breaking chiral invariance for k = ||. If 
we neglect thermal contributions to the momentum dependence of the fermion two-point function, or in 
fact for any truncation that assumes analyticity for the momentum dependence of the fermion two-point 
function, we will miss this important finite temperature contribution. The effect is numerically small in 
the presence of a large enough chiral symmetry breaking mass, but in the situation discussed above with 
k = 0, this effective mass-term takes care of the decoupling of the fermionic fluctuations, as we will argue 
in the following. 

Let us first briefly review the structure of the leading contribution to the fermion self-energy for large 
T in perturbation theory Q. The contributing diagram is given in figure 1. In the framework of the 
hard thermal loop approximation one then neglects the external momentum squared, p 2 , as well as the 
masses of the particles in the loop compared to the loop-momentum squared. After these simplifications 
the integrations may be performed and one obtains for the two-point function of a massless fermion 
(p = v ^2) 



(AfyHp) 



T 2 



1 



P0 In 1 P0 + P 



where m 2 -, is a quantity with the dimension of (mass) 2 , given by 



h 2 



2P \p - P 
P 2 V P 2 



Po + P 
Po-P 



(52) 



(53) 



The dispersion relation corresponding to ( |52|) is given by (we use u = (1,0,0,0) for the rest-frame of the 
heat-bath as is commonly done) 



= D{lo) = [1 + a(u, P)f (to 2 - P 2 ) + 2 [1 + a(uj, P)] b(w, P)uj + b(u, Pf 



(54) 
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This equation has four solutions ±u>± which can however only be given analytically in limiting cases. A 
general feature is that the solutions u = ±u>+ have the same (positive) ratio of helicity over chirality and 
are thus connected to the usual solution at vanishing temperature. The other solutions, u = ±w_ have a 
negative helicity over chirality ratio and constitute a collective mode termed the "plasmino". As P — > 0, 
both solutions approach each other and 

P 

u±(P) m ~^ P m T ± — . (55) 

In the opposite limit, for P S> %, one finds that uj^(P) exponentially approaches P with an also 
exponentially suppressed residue, whereas lu+(P) — ► P + rn^/P. The most important observation is 
however that the complicated momentum-dependence of the effective fermion propagator as given in ( |5"2| ) 
gives rise to a dispersion relation which resembles for small momenta P that of a massive particle with 
a mass oc T, while chiral symmetry is not broken by this behavior. It is also clear that the result given 
in ( |52| ) may not be reproduced using a derivative expansion as is routinely done in calculations using 
variants of the Wilsonian renormalization group. Since on the other hand in the framework of the TRG 
it is exactly this behavior that guarantees fermion-decoupling, we need to find a way to incorporate the 
above results in an approximation-scheme for the flow-equations. We now turn to this question. 

The canonical way to incorporate nontrivial momentum-dependence of n-point functions in the 
renormalization-group approach would be to expand the effective action in powers of the fields rather 
than in powers of derivatives. In order to obtain a closed system of flow-equations one then has to trun- 
cate the infinite tower of n-point functions by hand, either by setting the n-point functions with n larger 
than some n max to zero, or by using e.g. perturbative expressions for those. We could in principle follow 
such an approach, using scalar n-point functions derived from an effective action as given above and 
not imposing any restrictions on the momentum-dependence of the fermion two-point function. Since 
however the scale-dependence of the fermion self-energy is given by a flow-equation which involves a loop 
integration over a fermion propagator (it is basically given by the A-derivative of the diagram in figure 1 
where the two-point functions are replaced by the modified two-point functions as introduced in section 
2), we would end up with a rather complicated integro-differential equation. 

On the other hand we may also use the knowledge about the structure of the leading thermal correc- 
tions to the fermion propagator as discussed above. This structure is not changed by the presence of an 
infrared-cutoff scale A. All dependence of the leading correction on A is in the thermal mass-term ( |53|) , 
which becomes a function of A. Thus instead of (f3^), the thermal fermion two-point function for finite 
values of A reads to leading order in T, taking into account a possible chiral symmetry-breaking mass 
m xS B now, 



(A^-^A) = (l + a(po,P;A))jrf + b( Po ,P;A)^ - m xSB 
a( Po ,P;A) 



m 



T,A 



p2 



Po ln (PO + P 



b( Po ,P;A) 



rn 



P 



2P \ Po - P 



P 2\P 2 J \po-P 



(56) 



Correspondingly, the dispersion relation is obtained from (|54j) by replacing mj- by mr,A an d subtracting 
m xSB- The solutions have the same general features as discussed above. If we use this approximation 
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for the fermion propagator in the flow-equations for the effective potential and the Yukawa-coupling, we 
find the following result: 



Ad A U'( P ) = 



A 3 

-2N c \h 2 

7T 



N(u>+) N(u-) 



(57) 



and 



Ad\h 



A3_ 

"4vr 2 



u D(u ) 



+ 2 



D'(Q 



+,o) w 



-,o 



(58) 



Here uj n , and wo have been given in ( |36|) and (|38|) , whereas w± and u;± : o are the positive solutions of 

(59) 



U2 

= D(lo) = [1 + a(u>, A; A)] 2 (w 2 - A 2 ) + 2 [1 + a(u>, A; A)] 6(w, A; A)w + 6(w, A; A) 2 - 



at p ^ and p = respectively. Also primes on denote the derivative of D(u>) from (|59| ) with 

respect to w. We now obtain an additional flow-equation which governs the dependence of the thermal 
mass rriT.K (appearing in a and b above, see (56)) on A. This flow-equation reads 



Ad\m 



T,A 



-h: 



A 3 (n(lj ) 



4vr 2 



V 



+ 2 



D'(u +!0 ) D'(w_, ) 



(60) 



In this way we again have a closed system of differential equations which furthermore - in contrast to (|35| ) 
and (^), being obtained in lowest order derivative expansion - contains the most important corrections 
induced by thermal fluctuations. Even though the boundary value for the thermal mass term niy A for 
A — ► oo is zero, by (60) for any finite A > and h 2 one gets a nonvanishing effective thermal mass for the 
fermions even for p = 0. In this way, remembering the discussion above, fermion decoupling is guaranteed 
also at the second order phase-transition, since again the "energy-denominator" measured in units of A 
diverges as A — * 0. 

We could thus stop the discussion at this point and use equations (j57|) - (|^) to discuss the critical 
behavior of our model. However, we shall now show that following a proposal made in |30[| (in a different 
context), one may furthermore simplify the system of flow-equations without introducing large errors 
in interesting quantities. This further simplification will turn out to be very useful once we aim to go 
beyond the lowest order in the derivative expansion in the scalar sector. 

The key observation is that one may greatly simplify the effective fermion propagator without changing 
the most important features by writing 

(1 + o(po, P))j) + b(p , P)i + m xSB 



(1 + o(po, P)) V + 2(1 + a(po, P))b( Po ,P)po + b(po, P)< 
i> +m xSB + 0{m 2 T /P) 



m 



xSB 



p z 



m 



Isb + 0(m^P 2 ) 



(61) 



and neglecting terms suppressed by powers of P. Clearly, for m x sB = 0, the resulting expression preserves 
chiral invariance on the one hand, but displays the behavior of a massive propagator as far as the 
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0.05 



Figure 2: The dimensionless minimum kq as a function of t at the critical temperature for various values 



of x 



"(T=0) 
9(T=0) 



(x = 0.01, 0.2, 0.4, 0.6, 0.8, 1, 1.2, 1.4 from lowest to uppermost curve). 



dispersion relation is concerned on the other hand. If we use this modification, the flow-equations for the 
potential, the Yukawa-coupling, and the thermal mass read 



Ad A U'( P ) = 



W"( P ) + 2pU^(p)) ^1 + 3U"(p) 



2 



N c —h 



TT 



Ad A h 



A 3 



4vr 2 U'(0) - 2m 2 



T,A 



iV(cuo) | NjO^o) 



(62) 



and 



4vr 2 



jV(go) | jV(^, ) 



wo 



(63) 



with u^l = A 2 + 2m 2 TK + m 2 xSB , 2 fi 



A 2 + 2m| 1A . Since the fermions decouple from the universal critical 
behavior both in the simplified version of the flow-equations ((|62|) and (|63[)) as well as in the full HTL- 
improved versions (^)-(^), all universal results will be independent of the simplifications performed in 
order to arrive at (|62^). 

Before we turn to a discussion of the universal critical behavior in the next section, let us consider 
the modifications brought about by the emergence of a thermal mass for the fermions and the quality of 
the approximation (^). In order to check the universality of the critical behavior we may consider the 
model for different values of the Yukawa coupling h^ T=0 ^ at the respective critical temperature. If we 
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Figure 3: The value of the thermal mass m\ A divided by A as a function of t at the critical temperature 
for various values of x (x = 0.01,0.2,0.4,0.6,0.8, 1, 1.2, 1.4 from leftmost to rightmost curve). 

approach the critical temperature from below, the effective potential always has a nontrivial minimum 
at a finite value of n = Kq. If we consider the rescaled field and couplings according to ([4l"|), we expect 
the value of kq to be independent of A for small values of the flow-parameter. In figure 2 we display the 
dependence of the rescaled minimum of the potential on t = In A at the critical temperature for various 
initial values of the Yukawa coupling (unless stated otherwise, we always choose <?(t=o) =0.1 and vary 
x = ^?y =0 )/ff(T=o))- A scaling solution is reached for t % —7 and one nicely observes the independence 
of the final value ko,* of x. It is remarkable that the fermion-contribution is sizable down to t < 0, 
i.e. A < T. 

In order to discuss the importance of contributions neglected according to ( |61~D we display in figure 3 
the value of the thermal fermion-mass A divided by A (all dimensionful quantities are given in units 
of yjp(T=Q))- This is the relevant quantity in applications of (|6l]) in the framework of the TRG with a 
sharp cutoff, since the three-momenta inside the loop contributing to the flow of any quantity are P = A 
due to the J-function in the kernel ((|3l|) and correspondingly for the fermionic degrees of freedom). We 
plot in figure 3 m^ A /A again as a function of t for the same values of g^T=o) an< ^ x as use d in figure 2. 
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Figure 4: The thermal mass m| a-> in units of T 2 at the critical temperature as a function of x (solid 
line). For comparison the leading order result from ( |53[) is also given (dashed line). 



We note that the thermal mass remains very small for large t — a one-loop calculation for finite A yields 

1 1 



rriTA 



dxh 



A 
12 

R. 

2vr 2 



T 2 + 



l x 
4tt 2 

/i 2 



+ 



4vr 2 



T 



e x/T _ 1 e x/T + 1 

Tdilog(e A/T ) + Tdilog(e A/T + 1) + Aln(e A/T + 1) 



Te' A/T (A + T) (A/T large) 



(64) 



and rrij, A /A is smaller than 0.2 for the values of x considered here for all t^, — 2.5. On the other hand, as 
is obvious from figure 2, the main contribution from the fermionic degrees of freedom is at larger values 
of t. Thus we believe that the simplification of the fermionic dispersion relation made here introduces 
only very small errors for the nonuniversal results - as was pointed out above, universal quantities are 
completely unaffected^. 

Before turning to the phase-transition itself, consider finally the value of A — > 0. This is 

displayed in figure 4, where we have plotted m\ A ^ /T 2 at the critical temperature as a function of 
x again. For comparison the one-loop result (|53|) is displayed as dashed line. The renormalization of 
the Yukawa-coupling through thermal fluctuations yields a small but non-negligible deviation of this 
*Also note that neglected contributions in ( |6l| ) are formally of higher order in h 2 through (^). 
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Figure 5: The Yukawa-coupling hh—,o(T c ) as a function of x (solid line). For comparison h(r=o) 1S also 
given (dashed line). 

quantity from its lowest-order value. We point out that here the HTL-contribution (in the fermionic 



5 The universal critical behavior at the chiral phase-transition 

In this section we finally discuss the results for the chiral phase-transition as obtained from the quark- 
meson model after including the leading thermal effects in the fermionic sector as discussed above. 
We have pointed out that the derivative expansion fails for the fermion propagator and that fermion 
decoupling at a second order phase-transition is only guaranteed by the nonlocal fermion mass-term 
produced by the thermal fluctuations. After taking this effect into account, we expect the phase-transition 
to be in the universality class of the 0(4)-model. Critical exponents and amplitude-ratios should be 
independent of the number of colors or the Yukawa-coupling. Before we turn to the scalar sector and 
the universal quantities, let us briefly discuss two nonuniversal aspects: the critical temperature and the 
behavior of the Yukawa-coupling. 

In this section we introduce an additional effect, namely we go to first order in the derivative expansion 
in the scalar sector by taking into account a scalar wave-function renormalization. This quantity was 
already included in T\, (p8[). Whereas Z\ is in principle a function of the order-parameter field, we 



sector approximated by (|B1|)) is consistently resummed. 
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evaluate it at the minimum of the effective potential. Its i-dependence is given by 

dZ A 



dt 



-r]Z A 



(65) 



with the anomalous dimension r/. Having introduced Z\ we consider instead of the dimensionless quan- 
tities from (E3) renormalized dimensionless couplings according to 



Z^h 2 



n M = Zl n u {n) 



(66) 



We evaluate the anomalous dimension by considering the flow-equation for the 7r-two-point function 
and expanding it around q 2 = 0. Since the heat-bath introduces an additional four- vector u M , there are 
in principle two different wave-function renormalizations for (choosing as usually = (1,0,0,0)) q and 
qo. We do not make this distinction and consider only the effect of Z\ connected to q 2 , since this is the 
quantity which through ( |65|) yields the critical exponent ry (see below). For the anomalous dimension we 
find the following expression: 

A — — I F(oj a ,uj n ) + F{lo 7 ,,uj (7i 



+N C 



Utt 2 
h 2 



12tt 2 



[6uj^ - 3iv 2 - 9htu; 2 + 10/x 2 -] (-n(Au^) + Xu^h(Xu^,){h(Xuj^) - 1)) + 

[ufy + 3uj 2 Ht ~ 4 ^r] X 2 J 2 h h(Xio^){n(Xuj^) - l) (2h(Xuj^) - l) + 
^/xf,A 3 u^,n(Au;^) (rt(Aaty) - l) (6n(Au^) 2 - 6n(A^) + l) j (67) 



with 



F(w CT ,a; 7r ) 



5— -3 { n{Xuj a )(6ujl - &4 - 6 <4<4) 

+n(Aw„)(3<4 - - 10i^ 2 + 12u 3 a u>* + 15^^ - 9w ff a;£) 

-An(Aw w )(n(Au^) + l)(3^w 3 - 307*0/* + 10u/ 3 o/ 3 - QujIujI + 3w CT ^ - 7a/ ff c4; 
+A 2 n(Ai^)(n(Ac^) + l){2n(Xuj n ) + l)(w^ 2 - 2a£t4 + 



(68) 



2 

where /4 = /z£ 

renormalized quantities according to (^) should be used). In the limit where A is small compared to all 
other scales in the theory (this is the limit relevant for the fixed-point behavior) this rather complicated 
expression reduces to 



h K r fi 



^ = 1 + + Mx an d w ff as well as u> n are defined in (45) (where the 



V 



K T .( U {2) ) 2 4(4 - 3^1 - 3L0 2 a L0i + 44 



6vr 2 



+ 0(X) 



(69) 



and we again note the decoupling of the fermionic degrees of freedom in this limit. 

The flow-equations for the renormalized couplings introduced in (|6"6|) receive additional contributions 
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proportional the anomalous dimension. We thus have for example 



d t h 



-h + 



h 3 



(0) 



n(\u} r0 ) 



+ h(l) 



(-2 + rj)uy-> + (1 + rj)KrU. 



(2) 



A 

1^2 



+ 1k t u\, ' 



n(Xuj a ) 



+ 



n(Aw„ 



UJ-n 



and for the renormalized dimensionless minimum 



d t K, 



r,0 



-(1 + 77W0 + 



A 

4^2 



3 + 2/c r 



n(Xco a ) n(Xco n ) 



+ 3- 



to* 



+ 



JVc 2 n(Ao^) 

7T 2 Ofy 



(70) 



(71) 



Since the value of the anomalous dimension at the critical temperature will turn out to be very small, 
the modifications of the results displayed in figures 3 and 4 are tiny. 

The value of the unrenormalized Yukawa-coupling h 2 at T c as A —* is displayed as a function of 
x in figure 5. Also plotted is the (T = 0)-Yukawa-coupling and we note that it is only slightly modified 
through thermal fluctuations (for small values of the initial couplings). Here the introduction of a 
renormalized Yukawa-coupling however has important consequences: The value of rj is positive at the 
critical temperature. This leads to a vanishing renormalized Yukawa-coupling at T c , since in the infrared 
the flow-equation for h 2 reduces to 



d t h 2 



rjh 2 



(T 



T C ,A 



0) 



(72) 



Since 77 is very small, h 2 only vanishes very slowly (h 2 oc X v ). 

Let us now turn to the critical temperature. In leading order perturbation theory, one obtains T c as 
a function of x = h 2 T= ^ / g^T=o) as 



jnport 



4po(T = 0) 
1 + x 



(73) 



This dependence on x is basically reproduced by our resummed calculation. We display the critical 
temperature in units of a/ po(T = 0) as a function of x in figure 6. Note that in the limit x — > the result 
does not approach ([73]). This is of course due to the fact that this limit corresponds to a pure scalar 
theory where there are still modifications of T,P er [p4j. Indeed, the relative difference of our result with 
lowest order perturbation theory as displayed in figure 7 is small but nonvanishing for x — > 0. 

To end the discussion, let us report the results for the universal behavior. As the universal quantities 
of the 0(A^)-models are well studied for various values of N (for a review see [[HJ|) and the cases N = 1 
and N = 4 have been studied with the methods used here in pi] , 23, 24] we will only briefly discuss the 
main differences of our results with the ones that may be found in the literature. 

The universal critical behavior is encoded in the scaling form of the equation of state, which relates 
the order parameter p, the temperature T and the external field H = ■\j2pU\p) and may for example be 



written as [31] 
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Figure 6: The critical temperature in units of po (T = 0) as a function of : 
here). 
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(we use g(r=o) =0.1 



with the Widom scaling function f(x) and critical exponents 5 and 0. The Widom scaling function for the 
0(4)-model was calculated in the framework of the thermal renormalization group in [24], where however 
the anomalous dimension was neglected. It has also been computed from the exact renormalization 
group in the Matsubara- formalism pjjfl, from three-dimensional lattice simulations |p^|, and from mean- 



field theory as well as in second order e-expansion [33]. The TRG-results have been compared to those 
of the exact renormalization group in [24] and good agreement was found apart from the behavior for 
large x and T > T c . The disagreement in that region was due to differences in the value of the critical 
exponent 7, which governs the behavior of f(x) as x — > 00. In [24] it was argued that these differences 
should mostly be due to the neglected effects of the anomalous dimension. We are now in a position to 
check this claim through a calculation of critical exponents taking into account r\ according to fl67|). 

We have calculated the exponents (3 and 5 as introduced in (|74]), r\ as given in (|67|) , as well as 7 and 
v. The exponents 7 and v encode the behavior of the unrenormalized and the renormalized mass of the 
order parameter field as T — > T c according to 
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(75) 



whereas rj describes the behavior of the two-point function of the order-parameter field at the critical 
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Figure 7: The relative deviation of the critical temperature from the perturbative result, 5t c = 2 
as a function of x (g(r=Q) = 0.1). 
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temperature. Our results for these exponents are given in table 1, where we also give the results obtained 
from other methods for comparison)^. 

A useful consistency check is provided by the scaling relations which connect the critical exponents 



and leave only two of them independent. Thus for example one has [31] 

, d + 2-r] 

o = 

d-2 + rj 

7 = z/(2-r?) 



P 



(d - 2 + t?) 



(76) 



These scaling relations are fulfilled by our results to better than 0.5%. 

The results from |B5| are obtained from a fixed-dimension calculation in the three-dimensional 0(4)- 

model to seven loops whereas the results given in the last line are from three-dimensional Monte-Carlo 

simulations of the 0(4)-model [36|. We note that the results obtained in the present work, derived using 

the next-to-leading contributions in the derivative expansion in the scalar sector, are somewhat closer to 

the values found in lattice and three-dimensional perturbation theory-calculations than the results given 

^For results from the ERG see also where results that are essentially identical to those given in where obtained 
from a polynomial approximation to the effective potential. 
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0.429 


1.68 


0.85 


4.90 


0.017 


TRG + LO DE @ 


0.433 


1.73 


0.86 


5.0 


- 


ERG |2S] 


0.407 


1.548 


0.787 


4.80 


0.0344 


3d PT |H 


0.38 


1.44 


0.73 


4.82 


0.03 


3d MC [36| 


0.384 


1.48 


0.748 


4.85 


0.025 



Table 1: Critical Exponents for the 0(4) model. 



in [24|. However, the improvement for the exponents 0, 7 and v is perhaps surprisingly small for the 
present model. 

At this point it is interesting to compare the situation with the results from a study of the pure scalar 
0(l)-model. In j3?J a sharp-cutoff Wilsonian RG was used to obtain critical exponents of the three- 
dimensional 0(l)-model. We have explicitely verified that our flow-equation for the effective potential 
and the anomalous dimension, adapted to the O(l) scalar theory^ reproduce the expression given there 
in the scaling limit. The truncations made in |}7| are then identical to the ones made here and we 
would reproduce their values for the exponents. In that case, the exponents from the Wilsonian RG are 
in very good agreement with the results from other methods. This shows that in the 0(l)-model, the 
next-to-leading order derivative expansion as used here gives accurate results as one would expect given 
the smallness of the anomalous dimension^. The fact that the agreement in the present case is not as 
good might be due to the existence of a contribution in the present order of the derivative expansion 
which exists for N > 1 and has been neglected here, namely a term of the form Z\ip a ip a d 2 ifb^Pb (see also 
j24j). This question remains to be settled in future work. 



6 Summary and Conclusions 

In the present work we have extended the formulation of the thermal renormalization group '^M to 
theories involving fermionic degrees of freedom. Whereas the adaption of the full flow-equation (section 
2) is very much straightforward and does not involve any conceptually new points, it turns out that in 
performing approximations to the full functional differential equation one has to be considerably more 
careful for fermionic than for bosonic fields. After describing in section 3 the specific model for the Nf = 2 
chiral phase-transition that we studied, we have discussed at length questions concerning the decoupling 
of the fermionic fields at the phase-transition in section 4. 

The model under consideration consists of an 0(4)-symmetric scalar field coupled via Yukawa- 
couplings to 2N C chiral fermions. At low temperatures the 0(4)-symmetry is spontaneously broken 
down to an 0(3), where correspondingly 3 massless scalar fields ('pions') appear. The massive scalar 

*This is trivial for the flow-equation for the potential. The anomalous dimension has to be replace by the one defined 
from the Higgs since there are of course no Pions in the 0(l)-case. We have calculated the corresponding expression also in 
the present model. 

§ As a caveat, one should however notice that the scaling relations ( [7^ ) are violated by up to ~ 6% for the values given 
in @. 
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field is the would-be 'sigma', the fermionic fields acquire masses and the chiral symmetry thus is sponta- 
neously broken. This model is expected to be in the universality class of the three-dimensional 4-vector 
model with a second order phase-transition with well known universal behavior. The fermionic fields 
should not play a role for the infrared-behavior at the phase-transition and the effective potential should 
display scaling behavior for small field- values. 

We have then pointed out that the usual arguments about fermion decoupling at high temperatures do 
not apply in the TRG- approach in the presence of second- or weakly first order phase-transitions. The fact 
that the scalar distribution function dominates the fermionic one for small energies (or large temperatures) 
does not directly imply a dominance of the scalar contributions to the renormalized couplings (or the 
effective potential). The reason for this is the fact that it is not the distribution functions alone that 
enter the calculations, but the product of couplings times distribution functions. In the presence of a 
scaling solution for the effective potential, the couplings are also scale dependent in such a way that the 
product of scalar self couplings times scalar distribution functions has the same infrared behavior as the 
fermionic contribution, since in the latter the couplings do not exhibit scaling behavior. 

This is in contrast to the high-temperature limit, where neither the scalar nor the Yukawa-couplings 
scale. In this case it is indeed the relative behavior of the distribution functions that matters and fermion 
decoupling in the limit T — ► oo works by the usual arguments. 

One of the main points of the present work was then made in section 4 where we have discussed the 
question how universality in the presence of chiral fermions is nevertheless obtained and how this effect can 
be successfully implemented within the thermal renormalization group-approach. The key observation is 
that even in chirally symmetric theories thermal fluctuations modify the fermionic dispersion relations 
such that the corrected dispersion relation resembles that of a massive particle for small momenta [27]. 
This is the leading effect of finite temperature in the fermionic sector of the theory and is in this sense 
on equal footing with the scalar mass-correction oc T 2 which drives the phase-transition. Due to the 
chiral nature of the fermions, the thermal mass in the fermionic sector is however necessarily an effective 
mass resulting from a complicated momentum dependence of the two-point function. We have pointed 
out that this effect is not included if one applies a derivative expansion to the fermion two-point function 
as is generally done in the application of Wilsonian renormalization group-equations. While universal 
quantities are unaffected, the presence of the effective thermal fermion mass should modify nonuniversal 
quantities as obtained within the derivative expansion 2£j. 

We also discussed in section 4 how one may approximate the fermionic dispersion relations such that 
the resulting expressions are simple enough to keep the system of flow-equations manageable on the one 
hand and not loose the important effects on the other hand. Here we follow a proposal put forward in 
|30[| to rewrite the fermion propagator such that the thermal mass-term is effectively local, but chiral 
symmetry remains unbroken by it. We argue that this approximation is quantitatively reliable for small 
(T = 0)-couplings. 

After having clarified how fermion-decoupling works within the TRG we have finally discussed the 
results of a study of the chiral phase-transition in the framework of the quark-meson model in section 
5. Indeed, the phase-transition is second order with universal critical behavior as obtained from the 
pure scalar 0(4)-model in three dimensions. We find that for small enough (T = 0)-couplings the critical 
temperature is rather well described by perturbation theory. This observation was already made in 
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fl3| , 24, 38 1 for pure scalar 0(iV)-models. We improved upon previous results from the TRG on the 
critical exponents by including a scalar wave-function renormalization, i.e. going beyond leading order in 
the derivative expansion in the scalar sector. The critical exponents are displayed in table 1. We have 
also discussed the quality of our approximations as far as the critical exponents are concerned. It turns 
out that the inclusion of the anomalous dimension has a rather small effect on the exponents /3, 7 and v. 

In the present paper we have been concerned only with quantities which may also be obtained from 
Wilsonian renormalization group-calculations in the Matsubara-formalism. In principle, the formalism 



set up here is directly applicable to quantities like damping rates etc |25]. After we have clarified 
the necessary ingredients in the fermionic sector, we are now in a position to study out-of-equilibrium 
properties of the present model in linear response theory. 
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